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ABSTRACT 

We present distance estimates for eleven peculiar Virgo cluster spiral galaxies based on mea- 
surements of the stellar kinematics of their central 2 kpc. Stellar circular velocities were obtained 
using two-integral dynamical models. Distances were obtained by comparing, at each radius, 
the stellar circular velocities with synthetic Ha rotation curves derived from NIR Tully-Fishcr 
relations. The results show that most of our galaxies are located within 4 Mpc of the core 
of the cluster. Three of these galaxies, previously classified as "low rotator galaxies" or with 
"Truncated/Compact" Ha radial distributions, have stellar kinematics-based distances that are 
discrepant with Hi-based distances by at least 60%, and are likely to be located within the virial 
radius of the cluster. These discrepancies appear due to very truncated gas distributions plus 
non-circular gas motions or gas motions not in the plane of the stellar disk, perhaps as the result 
of gravitational interactions. Our results show that environmental effects can significantly re- 
duce the measured HI linewidths for some disturbed cluster galaxies, thus affecting the accurate 
determination of distances based on gas kinematics methods. 

Subject headings: galaxies: clusters: individual (Virgo) — galaxies: distances and redshifts — galaxies: 
fundamental parameters (classification, colors, luminosities, masses, radii, etc.) — galaxies: ISM — 
galaxies: stars — - galaxies: kinematics and dynamics galaxies: nuclei — galaxies: evolution — - galaxies: 
interactions — - galaxies: peculiar 
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1. Introduction 

For many years, the Tully-Fisher relation (Tully 
& Fisher, 1977) has been used to estimate dis- 
tances to the Virgo Cluster. These studies were 
based on optical or near infrared photometry 
(Gavazzi et al. 1999), together with gas kine- 
matics from HI linewidths (e.g., Yasuda et al. 
1997) or Ha rotation curves (Rubin et al. 1999). 
These techniques have been useful in the determi- 
nation of distances to individual galaxies within 
the cluster, allowing attempts to determine its 3- 
D structure (i.e. Yasuda et al. 1997, Gavazzi et 
al. 1999, and more recently Solanes et al. 2002). 

Recent results show that there is a population 
of galaxies in Virgo with HI line-widths narrower 
than those expected from the mean Tully-Fisher 
relation derived for the cluster. We refer to these 
as "Tully-Fisher deviant" galaxies. Either these 
galaxies are foreground to the Virgo cluster, or 
they have peculiarities which make their observed 
gas line-widths poor tracers of galaxy mass, and 
hence luminosity. 

The Rubin et al. (1999) study of Ha rotation 
curves of 89 Virgo disk galaxies identified a popu- 
lation ( ~ 6% of total sample) which appears to de- 
viate significantly from the Tully-Fisher relation, 
showing anomalously "low gas rotation" velocities 
at apparently every radius, which they called "low 
rotators galaxies". Solanes et al. (2002) studied 
the 3-D structure of the cluster using the Tully- 
Fisher relation based on HI linewidths and, most 
recently, 13 HI deficient spiral galaxies with appar- 
ently large 3-D barycentric distances located far 
away in front and behind the core of Virgo (San- 
chis et al. 2004) . Seven galaxies were identified as 
being in the foreground of the cluster core (D < 
13 Mpc), supposedly in areas where the ICM is 
too tenuous to be responsible for the stripping of 
the ISM gas from these galaxies. Of these galaxies 
only three had large projected distances (9 > 7°) 
, whereas the rest had small projected distance 
(9 < 4°). Some of the Solanes's small HI line- 
width galaxies are also Ha Rubin "low rotator" 
galaxies. 

In subsequent studies by the Sanchis et al. 
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group (Sanchis et al. 2002, Mamon et al. 2004, 
and Sanchis et al. 2004), the possibility was 
analyzed that these 13 Hi-deficient objects had 
crossed the core of the cluster, losing their gas 
content as a result of ICM-ISM stripping. They 
concluded these galaxies lie farther away than 
the maximum radius to which these galaxies can 
bounce out, making core-crossing improbable, and 
thus difficult for ICM-ISM stripping to be the 
cause of the observed Hi-deficiency. These au- 
thors suggested alternative scenarios including the 
possibility of erroneous distance estimations, HI- 
deficiency caused by tidal-interactions, heating of 
the gas by mergers, and even questioned the claims 
for Hi-deficiency. 

The application of the techniques discussed 
above to HI deficient galaxies suffers from the 
possibility that environmental effects can reduce 
the spatial extent of the atomic gas, cause non- 
circular motions or other disturbances, thus re- 
ducing the measured line-widths and therefore 
causing derived distances to be underestimated 
(Guhathakurta et al. 1988 and Teerikorpi et al. 
1992). Indeed these small line-widths could be 
caused either by truncation of the gas disks by 
ICM-ISM stripping, or because the gas velocities 
are intrinsically small at every radius (e.g. Rubin 
et al. 1999), which can happen if the gas motion 
is non-equilibrium or non-circular, or the gas ge- 
ometry is not understood. 

These conclusions raise obvious questions about 
these objects, such as: Could these "Tully-Fisher 
deviant" galaxies be the result of an underestima- 
tion of the distances produced by environmental 
effects, leading to an artificial distortion of the true 
cluster shape? 

We have carried out an extensive study of the 
kinematics and morphology of 13 peculiar Virgo 
Cluster galaxies (Cortes et al. 2006; Cortes et 
al. 2008 in prep). Six galaxies within our sample, 
which includes four tentative foreground galaxies, 
are objects with small HI line- widths and therefore 
"Tully-Fisher deviant" . In order to answer these 
questions, in this paper distances to these galaxies 
are estimated through an approach that relies on 
stellar kinematica3 and is therefore less affected by 



Unfortunately, two of the original sample galaxies; NGC 
4457 and NGC 4698, are unsuitable for deriving their dis- 
tances using stellar kinematics since NGC 4457 is almost 
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environmental effects. 

This paper is structured as follows. We describe 
in §2 the galaxy sample. In §3, the method for es- 
timating distances via stellar kinematics is intro- 
duced. The detailed application of this method to 
our data, the results on the distances to individ- 
ual galaxies of the sample, and the nature of the 
"Tully-Fisher deviant" galaxies is studied in §4. 
A discussion about the existence of the putative 
foreground group of galaxies in Virgo as well as of 
the causes for these galaxies Hi-deficiency is pre- 
sented in §5. Finally, suggested improvements to 
the technique, summary and concluding remarks 
are found in §8. 

2. The galaxy sample 

The sample consists of eleven peculiar Virgo 
cluster spiral galaxies, spanning a variety of op- 
tical morphologies (Table [1]). Our sample includes 
more early type than late type galaxies, since most 
of the strongly disturbed cluster galaxies are early 
types. Morphological selection was made using the 
R and Ha atlas of Virgo cluster galaxies of Koop- 
mann et al. (2001), whereas the kinematical selec- 
tion made use of the published Ha kinematics on 
89 Virgo cluster spirals by Rubin et al. (1999). 

The Virgo cluster has about 110 spiral and 
lenticular galaxies brighter than 0.1 L* Q (Bingelli 
et al. 1985, BST), so we have observed about 12% 
of the brightest Virgo cluster spiral galaxies. The 
sample galaxies are spread around the cluster core, 
spanning projected distances ranging from a few 
hundred kiloparsecs to 2.5 Mpc (Fig. [1]). 

While the sample selection is not uniform, it 
is designed to include bright Virgo spirals whose 
peculiarities are most poorly understood, and to 
include representatives of the different Ha types 
identified by Koopmann & Kenney (2004). In 
choosing sample galaxies within a given Ha type, 
we gave preference to those with kinematical pe- 
culiarities. 

Following the use of the star formation classes 
introduced by Koopmann & Kenney (2004), we 
can group our sample galaxies as follow 



face— on, so its rotation curve is uncertain, and NGC 4698 
has complex kinematics (see Bertola et al. 1999; Falcon- 
Barroso et al. 2006 ; Cortes et al. 2008 in prep) which 
makes its dynamical modelling challenging, and escapes the 
scope of this work. 



• Normal galaxies: NGC 4651. 

• Truncated/ Normal galaxies: NGC 4351, 
NGC 4569, NGC 4580, and NGC 4694. 

• Truncated/ Compact galaxies: NGC 4064, 
NGC 4424, and NGC 4606. 

• Anemic or Truncated/ Anemic galaxies: NGC 
4293, NGC 4429, and NGC 4450. 

Although our main goal is to study galaxy mor- 
phological evolution in the cluster, we have at least 
six galaxies in the sample with small HI linewidths 
(NGC 4064, NGC 4351, NGC 4424, NGC 4580, 
NGC 4606, and NGC 4694), including four galax- 
ies that Sanchis et al. (2004) suggest may belong 
to a foreground group. 

3. Distance estimations by using stellar 
kinematics 

The determination of distances to Hl-deficicnt 
galaxies using HI line-widths is vulnerable to 
the reduction in the spatial extension of the HI 
disks resulting from environmental effects, and 
by the possibility that gas velocities be intrinsi- 
cally smaller at every radius. On the other hand, 
the use of stellar kinematics for distance deter- 
mination ("stellar kinematics based distances", 
hereafter SKB distances) has the advantage that 
it is in principle unaffected by physical processes 
other than gravitational processes, whereas as we 
have seen, the kinematics of the gaseous content 
can be affected by a number of non-gravitational 
physical processes resulting from environmental 
interactions (e.g. ICM-ISM stripping). However, 
stellar kinematics are far more difficult to measure 
and analyze than gas velocities, on account of the 
faintness of the stellar component as well as the 
collisionless nature of the stars which makes the 
pressure supported component important. The 
latter implies that rotation velocities are no longer 
fully representative of the gradient of the poten- 
tial. Star velocities however, are as good tracers 
of the mass of the galaxy as gas velocities are, pro- 
vided one is able to properly handle their physical 
complexities. 

We must therefore remember that; a) it is im- 
portant to take into account the contribution of 
the stellar velocity dispersion a to the real circu- 
lar velocity, and b) we have to restrict the Tully- 
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Fisher analysis to a region within the galaxy where 
the surface brightness is high enough to allow re- 
liable measurement of stellar velocities. 

We have devised a distance estimation method 
based on stellar kinematics rather than gas veloc- 
ities. This method consists of the evaluation of 
the circular velocity from the stellar kinematics 
through the use of two-integral dynamical mod- 
els (see Binney, Davies & Illingworth 1990; van 
der Marel et al. 1990; Cinzano & van der Marel 
1994; Magorrian et al. 1998; Cretton et al. 1999) 
which are then used in a Tully-Fishcr relation. 
The method can be summarized as follow: 

1) The stellar kinematics are modeled via two- 
integral dynamical models, to yield the true circu- 
lar velocity V c . 

2) Because this circular velocity is valid only 
within the region that we have measured reliably 
it cannot be used in a traditional Tully-Fisher re- 
lation which requires sampling to large galacto- 
centric distances. To resolve this we construct 
synthetic ionized gas rotation curves for galaxies 
spanning a range of absolute magnitudes, using 
near infrared Tully-Fisher diagrams at each ra- 
dius, as in Rubin et al. (1985). These synthetic ro- 
tation curves represent the typical rotation curves 
as a function of absolute magnitude for "normal" 
galaxies. 

3) Finally, we compare these synthetic rotation 
curves with the circular velocities derived in 1). 
The best match gives the absolute magnitude, and 
therefore the distance modulus. 

In the following sections we provide a more de- 
tailed explanation of the estimation of the circular 
velocity and the construction of the synthetic ro- 
tation curves. 

3.1. Estimation of the circular velocity 
from stellar rotation curves 

The Multi-Gaussian Expansion MGE formal- 
ism (Emsellem et al. 1994; Cappellari 2002) pro- 
vides one of the simplest ways to build reliable 
two- integral dynamical models. The MGE formal- 
ism allows us to parametrize the surface bright- 
ness of a galaxy as a sum of gaussians with major 
axis widths Ci, axial ratios q[, and luminosities 
Li, from which we can find a unique luminosity 
density v (Cappellari 2002), assuming axisymme- 
try. This approach offers advantages over the stan- 



dard bulge + disk decomposition as a description 
of the light distribution. It provides an analyti- 
cal and efficient way to do the luminosity depro- 
jection, and it allows kinematical quantities such 
as the gravitational potential, velocity dispersion, 
and circular velocity to be estimated with a sim- 
ple integration, thus simplifying the calculations. 
These models can succesfully reproduce the ob- 
served stellar kinematics in (presumably) axisym- 
metric galaxies. However, we must keep in mind 
that issues such as dust obscuration, triaxiality, 
kinematically distinct features in the stellar disks, 
and the assumed lack of dark matter can easily 
make these models discordant with the observed 
kinematics. The bias introduced by these issues 
are discussed in §3.3. 

To estimate the circular velocity, two-integral 
self-consistent models under the MGE formal- 
ism were built for eleven galaxies of the sample. 
These models assume axisymmetry and <jr = a z . 
They probably do not reproduce properly the kine- 
matics, since the galaxies can exhibit important 
anisotropics in their velocity ellipsoid, but for the 
purposes of this work these models are acceptable 
since our aim is to estimate the circular velocity 
rather than fully reproduce the internal kinemat- 
ics. It is important to notice, that these models 
are purely stellar and do not include the contribu- 
tion from dark matter. In most of our galaxies our 
measurements are restricted to the inner parts (~ 
20% of the optical radius) making it impossible 
to constraint properly the dark matter halo pa- 
rameters. Thus, these SKB distances must in all 
the cases represent a lower limit to the distances. 
However, precisely because we are sampling the 
central regions where dark matter is the least im- 
portant its effect on our distances is probably mi- 
nor. The detailed modeling proceeded as follows: 

1) The MGE formalism was used to find the 
luminosity density v(R, z) from R-band flux- 
calibrated images (Koopmann et al. 2001), assum- 
ing an oblate system. The dust lanes were care- 
fully masked to minimize the effects of dust ob- 
scuration. For a given surface brightness S(x, y), 
the formalism finds the best MGE model (Fig [2} , 
which can be de-projected to yield a unique lu- 
minosity density. The inclination was assumed to 
be the minimum required for an oblate distribu- 
tion. In most cases this inclination was consistent 
with the inclination obtained by Koopmann et al. 
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(2001). 

2) If we consider the axisymmetric case, follow- 
ing Cappellari (2002), the MGE mass density can 
be written as 



P (R,z) = rJ2 



L, 



— ^ (o-jV^r) 3 ?.. 



exp 



R z + - 



i 

(i) 

where, T is the stellar mass-to-light ratio, N is the 
number of adopted Gaussian components, each 
with major axis width Ui, intrisic axial ratio qi, 
and luminosity Lj. 

Also, the MGE gravitational potential can be 
written as 



$(i?,z) = ,/^GT^^ f 1 Hi(T)dT, 

V 7T ^(TJo 



(2) 



where, G is the gravitional constant, and TLi is an 
dimensionless function given by 



Hi{T) 



exp 



2o? 



R 2 



1-£,T 2 



(3) 



with €i = 1 — qf. 

The circular velocity at each radius in the plane 
of the galaxies was calculated directly from the 
potential. The models have as input parameters 
the mass-to-light ratio T and the distance d, but 
the luminosity scales with the distance as L cx d 2 . 
This lead us to conclude that the circular velocity 
obeys 

V c 2 (R)=R^ = aF(R)- (4) 
Here a = Td, and T{R) is 



T[R) = 



G 



N 

E 



Li R 2 



/ Hi(T)T 2 dT 7 
Jo 



(5) 

which represents the radial variation of the circu- 
lar velocity, and doesn't depend on the distance. 
Thus, the models can be rescaled to any desired 
value of a. 

3) The ^ 2 = y/vl os ^^ LOS 



moment was cal- 



culated, for any point in the sky, following equa- 
tions 61-63 from Emsellem et al. (1994), after tak- 
ing into account the typo error noticed by Cappel- 
lari et al. (2006). To include the effects of seeing 
in the computation of /12, we weighted the data 



when appropriate by the seeing convolved surface 
brightness. This moment has the advantage that 
along the major axis it does not depend on the 
shape of the velocity ellipsoid^. 

It can be shown that [i 2 , scales with a, which is 
the quantity we wish to determine. Thus, we can 
choose a model with an initial guess for a (called 
ao, see table [3] for the value of czq used in each 
case ) , based on reasonable guesses for the galaxy 
distance and T, and rescale it to any mass-to-light 
ratio T and distance d as 



(6) 



4) Finally, the observed [12 was calculated us- 
ing the stellar kinematics data from Cortes et al. 
2008 (Fig. [3]). The model-derived /Z2(ao) moment 
and the observed 112 were compared, and a was 
calculated along the major axis of the galaxy as 



a = ao 



(7) 



yielding the mean a factor (e.g. Fig 2]). Er- 
rors were estimated as the standard deviation in 
a. Therefore, the sought-after circular velocity 
in equation (4) corresponds to that of the model 
which best matches the observed \X2 giving the 
best a. It is important to notice that here we 
have derived a, but not T and d independently. 
These should be disentangled from a. 

3.2. The Tully-Fisher relation for radii 
smaller than i?(Vmax) 

So far we have derived the best circular veloc- 
ity using two- integral dynamical models. In what 
follows we focus on making use of these veloci- 
ties to obtain the distances to our galaxies. In 
placing our circular velocities on a Tully-Fisher 
plot to derive the distance and additionally dis- 
entangle T from the a factor, we must remember 
that we have measured velocities only until R ~ 
0.4 i?20.5 (where -R20.5 is the radius at which the 
galaxy has a surface brightness in H-band of 20.5 
mag arcsec -2 ). This makes it impossible, as dis- 
cussed, to place our velocities on the usual Tully- 
Fisher relation, usually determined at ~i?20.5j m 



2 Unfortunately, along the minor axis it depends critically on 
the ratio j3 = (Jr/cTz. 
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order to obtain reliable distances. To overcome 
this problem, we use the approach of building syn- 
thetic rotation curves based on Ha gas velocities 
as in Rubin et al. (1985). These synthetic rota- 
tion curves represent the mean rotation curves for 
galaxies with different absolute magnitudes. 

The procedure to derive the synthetic rotation 
curves was the following: 

1) We constructed 77 c -band infrared Pseudo 
Tully-Fisher (PTF) relations for different normal- 
ized radii (0.05, 0.1, 0.2, 0.3, 0.6, 0.7, 0.8, 0.9 
R20.5), using the Ha rotation curves for Virgo clus- 
ter galaxies obtained by Rubin et al. 1999. (Fig. 
[5]). The PTFs were determined using the H-band 
magnitudes (Table ^ derived by Gavazzi et al. 
(1999), and Jarrett et al. (2003) for the galaxies 
of the Rubin's sample. The advantages of using 
infrared magnitudes rather than B magnitudes lie 
in the reduction of the dependence of the Tully- 
Fisher relation on Hubble type, in the reduction 
of the spread due to dust absorption and stellar 
population variations, and in the improved ability 
to trace the luminous matter (e.g. Pierce & Tully 
1998, and Gavazzi et al. 1996). We excluded from 
the T-F sample those galaxies that were identi- 
fied as belonging to possible background groups 
(d > 23 Mpc) by Solanes et al. (2002), whose 
could have introduced undesirable systematic er- 
rors. We finally ended up with 34 galaxies, all with 
Ha kinematics and infrared photometry, which we 
used to build the PTFs. We supposed that all 
these galaxies were located at a mean distance of 
16.8 Mpc, which corresponds to the distance to 
M 87 (Gavazzi et al. 1999; Nielsen & Tsvetanov 
2000; Tonry et al. 2001). The usual standard 
Tully-Fisher relation corresponds to a PTF mea- 
sured at ~i?20.5- 

2) We fitted a straight line to each of the above 
PTF relations corresponding to a given radius. 
Then, for every magnitude, we derived an interpo- 
lated gas velocity foe each radius from the family 
of fitted straight lines (Fig [6l top panel). 

3) The final curves for each magnitude (and 
thus mass) were constructed using the gas ve- 
locities obtained for each radius (Fig [HI bottom 
panel). These synthetic curves represent the ex- 
pected shapes of the rotation curves for galaxies 
with different luminosities (mass) belonging to the 
Virgo Cluster, and obeying the Tully-Fisher rela- 
tion. As we see in the graph, the synthetic rota- 



tion curves do not overlap, depending only on the 
absolute magnitude. 

The absolute magnitude corresponding to each 
synthetic curve was calculated by assuming that 
each curve correspond to a galaxy located at a dis- 
tance Dy = 16.8 Mpc. Finally, the absolute mag- 
nitudes M of our sample galaxies could be found 
by locating the synthetic rotation curve that best 
matches the stellar circular velocity curve for each 
sample galaxy. The latter can be determined by 
minimizing the \ 2 (Figure [7]) between the stellar 
circular velocity and the synthetic rotation curves. 
The distance modulus follows from the observed 
apparent magnitude m. Errors were estimated us- 
ing the variation A% 2 , as 3-cr errors (A^ 2 =9). 

Having derived the absolute magnitude by 
matching the best model circular velocity and the 
synthetic rotation curves, the distance d is simply 
derived as 



Moreover, after deriving d we can obtain the mass- 
to-light ratio as T = a/d. 

The comparison between gas-based synthetic 
rotation curves and stellar circular velocities is 
justified, if we exclude the galaxies with pecu- 
liar rotation curves such as the low rotators, since 
gas-based synthetic rotation curves represent the 
typical rotation curve, and thus the gradient of 
the potential. It is important to notice however 
that discrepancies between gas-based synthetic ro- 
tation curves, and stellar circular velocities could 
be found in the inner parts of the galaxies (r < 0.1 
-^20.5)1 where the gas velocity dispersion or non- 
circular gas motions could be important enough to 
affect the determination of the gas-based synthetic 
rotation curves, and peculiar stellar structures as 
circumnuclear disk can be found. Thus we restrict 
our comparision between stellar circular velocity 
and gas-based synthetic rotation curves to radii 
greater than 0.1 i?20.5- 

3.3. Biases in the stellar kinematics-based 
distance estimations 

Our SKB distances are not free of systematic 
effects. We recall that these effects are: 

1) The presence of dust in our galaxies can 
reduce the amount of luminous mass that con- 
tributes to the potential, so we could underesti- 
mate the expected circular velocity. This problem 
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arises preferently in galaxies with prominent dust 
lanes. 

2) Triaxiality and non-axisymmetric structures 
are ignored in the pursuit of oblate models. This 
in turn introduces either an overestimate or under- 
estimate of the mass by altering the amplitude and 
shape of the rotation curve. Galaxies with strong 
bars could be seriously affected. Moreover, triax- 
ial features are very difficult to model and their 
inclusion is beyond the scope of this work. 

3) Kinematically distinct components are also 
important. Some galaxies have circumnuclear stel- 
lar disks (e.g. NGC 4429) or counter-rotating 
disks which can be easily detected in the veloc- 
ity fields but are hard to detect photometrically. 
The latter makes the inclusion of such components 
in the models difficult making our models discor- 
dant with the observations in regions dominated 
by such features. Thus, we must exclude in our 
comparison the regions where these features are 
important. 

4) The dark matter halo is an important com- 
ponent of the potential of the galaxy especially in 
the outer parts. The dark matter contribution can 
be represented by a logarithmic potential, 

<^M = f Mr 2 +r 2 ), (9) 

which is the simplest potential producing a flat 
circular velocity vq at larger radii (r >> r*o). In 
figure [SJ we show the effect of the dark matter 
halo in two sample galaxies; NGC 4450 and NGC 
4569. The parameters tq and vq were chosen in 
order to match the observed Ha rotation curve 
from Rubin et al. (1999). This figure shows that 
the dark matter halo is indeed important in the 
outer parts, but our observations are limited to 
the inner 30" , so we are unable to constrain the 
dark matter halo parameters properly, implying 
that we are underestimating the circular velocity. 

Most of these problems imply an underestima- 
tion of the circular velocity. For example, dust 
and dark matter have the similar effect of hiding 
the amount of mass that contributes to the poten- 
tial. This translates to an underestimate of the 
luminosity and therefore the model masses. Thus, 
we should consider these SKB distances as lower 
limits to the real distance; galaxies will tend to be 
somewhat further away than estimated here. 



4. Distances to the peculiar Virgo cluster 
galaxies 

4.1. The actual fittings 

Now that we have discussed in detail the 
methodology, we focus on the actual fittings and 
their errors. The MGE fits to the light distribu- 
tion are presented in Fig. [5] We see that most 
galaxies are reasonably well fit by the MGE algo- 
rithm. Differences in surface brightness between 
the galaxy images and the MGE model are less 
than 10%, which translates to velocity uncertain- 
ties of the order of 5-10%. Galaxies such as NGC 
4064 and NGC 4450 present bigger uncertainties 
resulting from the existence of a central bar that 
obviously cannot be accounted by an axisymmet- 
ric decomposition. 

The a/ao ratio salong the major axis are dis- 
played in Fig [U In most galaxies, the values are 
roughly constant along the major axis outside the 
central 5", except in NGC 4293, NGC 4424, NGC 
4429, and NGC 4450. In NGC 4429 there is a 
prominent circumnuclear disk which cannot be re- 
produced with these simple models. The best 
a/ao ratios, therefore a, were determined (Table 
[3]) as the average along the major axis outside the 
central 5" , except for NGC 4429 where we use dis- 
tances greater than 12" in order to exclude the 
circumnuclear disk. The average a/ao ratio is dis- 
played as a solid line in figured] (left panels), while 
the errors appear displayed as dashed lines. The 
comparisons between the observed and best model 
H2 moments along the major axis for the galaxies 
are displayed in FigH] (right panels) On these plots, 
the solid line represents the best model defined by 
the mean a factor, the shaded region represents 
the range of uncertainty given by the standard de- 
viation in a. Along the major axis, the model and 
the observed /12 hi most cases agree reasonably. 

To compare the synthetic rotation curves 
with the circular velocity curves, the latter were 
rescaled in radius in terms of -R20.5 in the H-band 
(Gavazzi et al. 1999). The exceptions were NGC 
4064 and NGC 4424 whose values were determined 
from H-band surface brightness profiles derived by 
Jarrett et al. (2003), since Gavazzi et al. (1999) 
do not have measurements for these galaxies. This 
is shown in Fig. [7] The model-based circular ve- 
locities extend further than the stellar kinematic 
data. This is justified since the circular velocities 
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are derived from self-consistent models based on 
optical images which extend into the outer galax- 
ies. 

The comparison between the stellar circular ve- 
locities and the synthetic rotation curves were per- 
formed by monitoring the x 2 between the syn- 
thetic rotation curves and the stellar circular ve- 
locities, within the range of 0.1 - 0.8 i?20.5- This 
range was chosen in order to avoid large variations 
in velocity existing within the inner 0.05 -R20.5 (e.g. 
NGC 4429, NGC 4450, and NGC 4569), and to 
minimize the effect of not including a dark mat- 
ter halo, which can be significant if we compare 
only the outer parts (i.e., using a range of radii of 
0.6 - 0.8 ^20.5)- Comparison between the present 
circular velocities with those derived by Rubin in 
Ha show that in most cases Rubin's values are 
of the order of or lower than our circular veloc- 
ities. The exceptions are NGC 4450, and NGC 
4569. In both cases, this implies that the galaxies 
are more massive than predicted by our models. 
In the case of NGC 4450, the MGE fit (Fig. [5} 
exhibits a big discrepancy in P.A on one of the 
isophotes, indicating the presence of an important 
non-axisymmctric structure which might account 
for the missing mass. In NGC 4569, dust lanes are 
prominent until at least 5 kpc through the disk, 
which prevented us from including all the lumi- 
nous mass. 

4.2. The SKB distances 

Once the absolute magnitudes M# for each 
galaxy is determined as is described in section §3.2, 
the line-of-sight distance can be found from the 
distance modulus. T can also be derived from d 
and a (Table [4]). To test whether the results on 
T are physically meaningful, we plot in Fig [9] the 
relation between the R-band luminosity, and the 
R-band mass-to-light ratio. The best fit straight 
line in logL, log Y space is 

log(T/T Q ) = -3.2±0.6 + (0.35±0.05)log(L/L o ) 

(10) 

This fit is roughly consistent within the errors 
with the fundamental plane correlation T cx L 2 
predicted from the virial theorem (e.g., Faber et 
al. 1987; Bender, Burstein, & Faber 1992). This 
consistency lends credibility to the modelling pre- 
sented here. The observed variations in T are 
mostly due to the variation in L. However, there 



are galaxies that present very low T. For example, 
NGC 4694 has a value of T of 1.1 ± 0.3 T Q , which 
is a factor ~ 2 lower than the expected. This low 
value and its B — V color (~ 0.62), suggest a recent 
starburst. 

Figure [10] which displays the galaxy distance 
distribution shows that our sample galaxies are 
located between 10 to 27 Mpc from us, with the 
majority of them between 13 and 20 Mpc. Two 
galaxies, NGC 4450 and NGC 4569, have model 
locations in front of the Virgo main core at about 
10-12 Mpc from us, although the models under- 
estimate the mass for these galaxies, so it is likely 
that these galaxies are not foreground (see notes in 
table ED . On the other hand, NGC 4651 is located 
at 27 Mpc and is clearly a background galaxy. 

It is interesting to compare our stellar kinematics- 
based distance estimates with Solanes's Hi-based 
distance estimates. In figure [HI we show the cor- 
relation between Solanes's distances and the ratio 
between the SKB distances and Solanes's, with 3- 
o~ error bars. ^From inspection of this figure, we 
conclude that NGC 4293, NGC 4450, NGC 4569, 
NGC 4580, and NGC 4651 have SKB distances 
that agree within the errors with those of Solanes. 
But others such as NGC 4064, NGC 4351, NGC 
4424, and NGC 4606 have significant discrepancies 
worthy of further analysis. 

4.2.1. The cases of NGC 4064, NGC 4424 and 
NGC 4694: Rubin's low rotators 

NGC 4064 and NGC 4424 are particularly inter- 
esting because both galaxies have been regarded 
as possible foreground galaxies. Moreover, both 
galaxies, together with NGC 4694, belong to the 
class of Rubin's "low rotator" galaxies, although 
for the case of NGC 4064 Rubin's Ha velocities are 
much smaller than our ionized gas velocity mea- 
surements (Cortes et al. 2006). Solanes's distance 
estimations indicate a distance of 9.8 ± 0.4 Mpc 
for NGC 4064, and 4.1 ± 0.2 Mpc for NGC 4424^1 
in stark contrast with our estimations of 18.0 ±Jg 
Mpc, and 15.2±1.9 Mpc respectively. The dis- 
crepancies are of a factor 1.8 for NGC 4064 and 
3.7 for NGC 4424, which are statistically signifi- 
cant. In the case of the other "low rotator" , NGC 
4694, we determine a distance of 13.4 ±\'q Mpc. 
Clearly, NGC 4064 and NGC 4424 belong to the 
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Virgo cluster. NGC 4694 might be affected by 
the dark matter bias in our models, so its dis- 
tance should be considered a lower limit. All these 
galaxies present truncated Ha disks (Koopmann & 
Kenney 2004), and their gas kinematics are pecu- 
liar (see Cortes et al. 2006), indicating that they 
are heavily disturbed by the environment. They 
all show evidence for strong gravitational interac- 
tions or mergers, which may account for their low 
gas velocities. This is cautionary of the risk of 
using gas velocities as distance estimators for HI 
deficient and/or disturbed galaxies. We conclude 
that the "low rotators" found by Rubin do not ex- 
hibit a low gas rotation velocity because they are 
intrinsically faint and thus closer to us, but rather 
because they have been heavily affected by grav- 
itational interactions perhaps in addition to ram 
pressure stripping. 

4.2.2. NGC 4351 and NGC 4606 

Consider now the cases of NGC 4351, and NGC 
4606. Both galaxies belong to Rubin's sample 
but they were not classified as "low rotators". 
Solanes's estimation gives distances of 11.2 ± 1.3 
Mpc, and 12.7 ± 0.3 Mpc respectively which again 
differ from our results of 20.3 ±?;| Mpc, and 19.9 
±H Mpc. In NGC 4351 the discrepancy, could 
be due to truncation of the HI disk via ICM-ISM 
stripping, a very likely effect on this galaxy con- 
sidering its high relative velocity with respect to 
the cluster mean (~ 1200 km s" 1 ). NGC 4606 
has a star formation class "Truncated/Compact", 
same as NGC 4064, and NGC 4424, and as these 
galaxies it has features that suggest gravitational 
interactions (Cortes et al. 2006), so again its small 
line- width could be due to the action of environ- 
mental effects. 

4.2.3. NGC 4569 

This galaxy has a SKB estimated distance of 
9.9 ±0.2 Mpc, still in front of Virgo. However, 
the SKB distance for NGC 4569 may be an un- 
derestimate, because of its strong dust lanes and 
possibly large influence of its dark matter halo 
(Fig [H]). If we include the dark matter halo in 
order to match the Rubin's Ha rotation curve we 
find a barycentric distance of 3.8 Mpc, locating 
this galaxy inside the maximum rebound radius. 
The SKB distance is slightly bigger than Solanes's 
Hl-based distances, although there are different 



HI linewidths reported in the literature for NGC 
4569, and one gets a different distance depending 
on which linewidths one uses. The HI linewidth 
in figure comes from the Arecibo General Cata- 
logue, a private database mantained by Riccardo 
Giovanelli and Martha P. Haynes, partially pub- 
lished by Solanes et al. (2002). This is not neces- 
sarly consistent with the homogenized Hi-distance 
of Solanes et al. (2002), which is an average be- 
tween different Tully-Fisher (TF) catalogues. If 
fact, the Hi-distance in the TF catalogues used by 
Solanes et al. can be as low as 7.9 Mpc (Kraan- 
Korteweg et al. 1988) and as high as 16.1 Mpc 
(Ekholm et al. 2000). 

4.3. SKB distances and star formation 
classes 

As it was shown in §4.2, in the cases of NGC 
4293, NGC 4450, NGC 4569, NGC 4580, and NGC 
4651, Hl-based distances agree with our SKB dis- 
tances. These galaxies have star formation classes 
(Koopmann & Kenney 2004) ranging from "Trun- 
cated/Anemic" for NGC 4293, and NGC 4450, 
to "Truncated/Normal" for NGC 4569 and NGC 
4580, and to "Normal" in the case of NGC 4651. 
In contrast for the "Truncated/Compact" galaxies 
the SKB distances don't agree with Hl-based dis- 
tances. This suggests that "Truncated/Compact" 
galaxies are not only very truncated in their gas 
content, but that they also have non-circular gas 
kinematics or gas motions which are not in the 
plane of the stellar disk. 

Here we comment on the possible existence of 
several galaxies in the foreground of Virgo at ^10 
Mpc (Solanes et al. 2002), which were tentatively 
identified as a group by Sanchis et al. (2002). 
Four of our sample galaxies are considered fore- 
ground by Solanes. All of them except for NGC 
4569 are "Truncated/Compact" (i.e., NGC 4064, 
NGC 4424, and NGC 4606) with SKB distances 
in the range 16-20 Mpc. NGC 4569 has a stellar- 
kinematics estimated distance of 9.9 ±0.2 Mpc, 
still in front of Virgo. However, this galaxy shows 
strong evidence for ongoing ICM-ISM stripping 
(Vollmer et al. 2004; Kenney et al. 2008 in prep), 
and its angular distance (1.7°) suggests that it lies 
within a region where the ICM medium is dense 
enough for efficient stripping of its gas. It is im- 
portant to recall that our SKB distances should 
be considered as a lower limit due to the biases. 
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We should expect that the real distances are larger 
than our estimates, so the discrepancies between 
Solanes' and our determinations should increase if 
we correct for the biases. 

We conclude, on the basis of our overlapping 
sample (we have analyzed 60% (4/7) of their fore- 
ground candidates), that there is no compelling 
evidence for the existence of multiple foreground 
galaxies at D~10 Mpc, as suggested by Solanes et 
al. (2002). 

5. Barycentric distances and the origin of 
Hl-deflcient galaxies in the outskirts of 
the cluster 

In order to determine how far these galaxies 
are with respect to the core of the Virgo Clus- 
ter, we calculated their barycentric distances r to 
M87, assumed to be at a distance of £V=16.8 Mpc 
(Gavazzi et al. 1999; Nielsen & Tsvetanov 2000; 
Tonry et al. 2001). These distances can be derived 
from 

r 2 = D v +d 2 - 2D v dcos6 (11) 

where d is the line-of-sight distance of the galaxies, 
and 6 the angular distance to M87. The results can 
be found in Table [4] and are shown in Fig [101 All 
galaxies are plotted on a P-V diagram. The solid 
line represents the "virial radius" (rv = 1-65 h^h 

Mpc, where h 2 / 3 = if /(66.7 km s _1 Mpc -1 )) de- 
rived by Mamon et al. (2004), and the dashed 
line represents the largest "maximum rebound ra- 
dius" {r re b = 1-2.5 ry) of all the sample galaxies. 
The latter corresponds to the maximum radius to 
which particles that cross the core can bounce to. 
In their calculations Mamon et al. (2004) con- 
sidered only one single halo, the M 87 subcluster, 
which is the most massive substructure in Virgo. 

We see in Fig [10] that most of our galaxies are 
within the largest maximum rebound radius, and 
such galaxies could have crossed the core of the 
cluster and suffered ICM-ISM stripping, which re- 
ducing the size of their gas disks. While many of 
the Hi-deficient spirals in the outer cluster prob- 
ably lost most of their gas during a core pas- 
sage, others have apparently experienced ICM- 
ISM stripping in the cluster outskirts. The Virgo 
spiral NGC 4522 is a clear example of outer cluster 
ICM-ISM stripping (Kenney et al. 2004, Crowl & 
Kenney 2006), which is thought to be caused by 
ram pressure which is enhanced due to an ICM 



which is dynamic and lumpy rather than static 
and smooth. 

Others likely owe their Hi-deficiency to gravita- 
tional interactions, such as NGC 4064 (Cortes et 
al. 2006) and NGC 4293, which exhibit kinemati- 
cal peculiarities and disturbed dust distributions, 
both signatures of gravitational interactions. 

NGC 4450, NGC 4569, NGC 4580, and NGC 
4651 are galaxies further out from the maximum 
rebound radius, according the SKB distances. The 
first three galaxies are HI deficient. NGC 4450 is 
an anemic galaxy with no evidence of recent ICM 
stripping. Its maximum circular velocity as deter- 
mined by us is lower than the maximum Ha veloc- 
ity, indicating that our model contains less mass 
than needed for reproducing the observed Ha ve- 
locities, probably due to either non-axisymmetric 
features which were not included or a significant 
dark matter halo. If wc include a dark matter 
halo in order to match the Ha rotation curve (Fig 
[8]), it is likely located inside the maximum re- 
bound radius (d ba r y = 3.3 Mpc). NGC 4580 has 
a "Truncated/Normal" star formation class, so it 
is likely to have been affected by ICM-ISM strip- 
ping. However it couldn't have crossed the core 
according to the simple cluster model, and the age 
of stellar population in outer disk reveals that it 
must have been stripped outside the cluster core 
(Crowl & Kenney 2008). This galaxy is closer (in 
the plane of the sky) to M49 than M87, so per- 
haps it is bound to the M49 subcluster rather than 
the M87 subcluster, and affected by the M49 en- 
vironment. The M49 subcluster is 4 times less 
massive than the M87 subcluster, and has a less 
dense ICM, so it is unlikely to strip the gas from a 
galaxy if the ICM is static and smooth, but could 
if the ICM were dynamic and lumpy. In general 
substructures in the cluster can be important in 
several ways: by stirring up the ICM through sub- 
cluster mergers and locally increasing the strength 
of ram pressure, by changing the value of the max- 
imum rebound radius, or by the pre-processing of 
galaxies in infalling groups. 

The case of NGC 4651 is clearer: this galaxy 
is very far from the core of Virgo (+10.3 i^'l 
Mpc), and its star formation class is "Normal", 
indicating that it has not been affected by the 
ICM medium, so it probably correspond to a back- 
ground galaxy rather than a Virgo member. 
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Finally, the case of NGC 4569 is the most com- 
plicated. The SKB distance places this galaxy out- 
side of the maximum rebound radius (r = -7.0 
±0.2 Mpc), and in front of the cluster. Its high HI 
deficiency, truncated HI and Ha distribution, and 
extraplanar gas arm (Vollmer etal 2004; Kenney 
etal 2008) are clear evidence for recent ICM-ISM 
stripping. Simulations (Vollmer et al. 2004) sug- 
gest that peak pressure occurred 300 Myrs ago. 
Its extreme line-of-sight velocity (-1200 km s _1 
relative to the cluster mean) is consistent with a 
galaxy after peak pressure located on the near side 
of the core, but it is hard to understand how it 
could be as much as 6 Mpc closer than the core. 
The prominent dust lanes are an issue here, bias- 
ing the distance determination, so it might simply 
happen that this galaxy actually lies inside the 
maximum rebound radius. Also, if we test the in- 
clusion of a dark matter halo in order to match 
the Ha rotation curve (Fig [5]), the corrected dis- 
tance is about 13 Mpc, putting this galaxy inside 
the maximum rebound radius. 

6. Conclusions 

In this paper, we have shown that stellar kine- 
matics can be reliably used for determining dis- 
tances to spiral and peculiar galaxies. This tech- 
nique presents a considerable advantage with re- 
spect to gas kinematics based techniques since, 
unlike the latter, it is in principle unaffected by 
physical processes other than gravity. 

This technique can be considerably improved 
provided we are able to overcome the biases de- 
scribed in §3.3: 

• Overcoming dust obscuration by using near- 
infrared imaging. 

• Extending observations of galaxies further 
out. This has the advantage of going be- 
yond bars, and bulges, and getting a better 
constraint on the dark matter halo. 

• Improving the reliability of synthetic rotation 
curves by increasing the number of sample 
galaxies used to build these curves, and ex- 
tend the sample to non-cluster undisturbed 
galaxies. 

Using stellar kinematics techniques, we have 
obtained SKB distances to 11 peculiar Virgo Clus- 



ter galaxies. Our results can be summarized as 
follows: 

1. - Of the four galaxies in our sample with HI 
line-of-sight distances estimated by Solanes et al. 
(2002) to be significantly foreground to Virgo, 
three galaxies have SKB distances within 4 Mpc 
of the cluster core, only one (NGC 4569) has an 
SKB distance (within 3<r) consistent with the HI 
LOS distance. 

2. - Galaxies with star- formation classes other 
than "Truncated/Compact" have SKB distances 
that agree with the Hi-based distances within 
3<7 errors. On the other hand, the "Trun- 
cated/Compact" galaxies NGC 4064 and NGC 
4424 exhibit large discrepancies between Hi-based 
and SKB distances. SKB predicts distances 2 to 4 
times farther away than Hi-based distances. Thus, 
environmental effects can indeed affect both their 
gas content and kinematics, rendering Hi-based 
distances unreliable. 

3. - SKB distance determinations show that "low 
rotator" galaxies belong to the cluster, strongly 
suggesting that environmental effects are the cause 
of the their low gas velocity. Galaxies with highly 
truncated gas disks and other peculiarities such 
as disturbed stellar disks and non-circular or 
non-planar gas motions give gas-based distances 
smaller than SKB distances. These galaxies were 
perhaps subject to gravitational interactions plus 
ICM-ISM stripping. On the other hand, galax- 
ies with normal or truncated gas disks, and nor- 
mal stellar disks, presumably affected only by 
ICM-ISM stripping, have SKB distances consis- 
tent within the errors. with gas based distances. 

4. - Most of the galaxies in our sample are within 
the bounce-back zone (e.g. distances within 4 
Mpc of M87), indicating they could have crossed 
the core of the cluster. If this is the case, ICM- 
ISM stripping could explain the existence of Hi- 
deficient galaxies in the outskirts of the cluster, 
although there are some galaxies that were clearly 
stripped outside of the core. This is probably due 
to the existence of a dynamic lumpy ICM stirred 
up by sub-cluster mergers, or gravitational inter- 
actions plus associated gas stripping. 
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13 h 00 m 12 h 45 m lS h 30 m 12 h 15 m lS^O" 1 
Right Ascension (JSOOO) 

Fig. 1. — Sample galaxies and their location in the Virgo cluster. Contour map represent the ROSAT 
X-ray emission in the cluster (Bohringer et al. 1994). Sample galaxies are represented by a black cross and 
its NGC name. Solanes's foreground galaxies are represented by a circle, whereas background galaxies are 
represented by a diamond. Four sample galaxies overlap with Solane's foreground galaxies, and one with 
Solanes's background galaxies. Black circle represent the virial radius (~ 5.6°) around M 87. 
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Fig. 2. — MGE fits to sample galaxy R-band images. MGE model contours are overplotted on R-band 
contours. 
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Fig. 3. — Stellar kinematics along the major axes of the sample galaxies. From left to right: the stellar 
velocity, the velocity dispersion, and the ^2 moment. 
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Fig. 3. — Continued. 
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Fig. 3. — Continued. 
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Fig. 4. — a/cto ratios (left panels) and /12 moments (right panels) along the major axes. The best fit model 
is represented by thick solid lines, and the model uncertainties are indicated by dashed lines (left) or shaded 
areas (right). 
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Fig. 4. — Continued. 
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Fig. 4. — Continued. 
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Fig. 5. — Pseudo Tully-Fisher relations for different radii. Solid circles represent Rubin's Ha velocities for 
Rubin's galaxy sample. The lines represents linear fits to the Pseudo Tully-Fisher relations. 
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Fig. 6. — Pscudo Tully-Fishcr relations and synthetic rotation curves. Top Panel: Lines fitted to the 
Pseudo Tully-Fisher relation for different radii (0.05, 0.1, 0.2, 0.4, 0.8 i?20.5)- From the interception between 
those lines and the line for constant magnitude, we find the corresponding rotational velocity for each 
radius. Bottom Panel: Synthetic rotation curves derived from Pseudo Tully-Fisher relations labelled with 
corresponding values of absolute magnitudes Mjj. 
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Fig. 7. — MGE stellar circular velocities and synthetic rotation curves. Left panels; stellar circular velocities. 
Circular velocities (thick solid line with shaded area representing uncertainty) roughly match with one 
synthetic rotation curve (solid lines). Ha gas rotation curves are represented by open symbols. HI maximum 
velocities (Solanes et al. 2002) are represented by the dashed line. The numbers to the right of the synthetic 
rotation curves represent the absolute H magnitudes for a galaxy having such a rotation curve at the distance 
of the Virgo cluster. Right panels; Ax 2 between the MGE circular velocity and synthetic rotation curves, 
dashed lines represent the 3-ct confidence level. 
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Fig. 7. — Continued. 
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Fig. 8. — Effect of the dark matter halo in two sample galaxies. Solid lines represent the circular velocities for 
different dark matter halo parameters. Open symbols represent the Rubin etal (1999) Ha rotation velocities, 
and shaded areas represent the uncertainties in the pure stellar circular velocities. Our stellar kinematics 
data is limited to the inner 30" , so we cannot constrain the dark matter halo parameters without using the 
ionized gas rotation curves, but even including them the dark matter halo parameters are too uncertain. 
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Fig. 9. — Correlation between the R-band stellar mass-to-light ratio T and the luminosity L. Solid line 
represents the fit of a straight line, which roughly follows the fundamental plane relation T oc L - 2 . 
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Fig. 10. — Position- Velocity diagrams for sample galaxies in the Virgo cluster. Top panel: LOS distances 
versus systemic velocities. Dashed line represents the systemic velocity of M87. Cross represents the location 
of the core of the Virgo cluster (M87). Bottom panel: Barycentric distances versus projected relative velocity 
with respect to M87. Solid line represents the virial radius, and dash line represents the maximum rebound 
radius. Arrows represent the expected location of NGC 4450 and NGC 4569 if a dark matter halo is included. 
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Fig. 11. — Ratio of SKB to HI based distances versus HI based distance. Galaxies with small LOS HI based 
distances tend to have SKB distances at least 50% higher. This shows that HI based distances can be biased 
by environmental effects. 
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Table 1 
Galaxy Sample Properties 

















Inc. 


P.A 




DM87 


Name 


R.A. (J2000) 


Docl. (J200) 


RSA/BST 


RC3 


SFC 




(deg) 


(deg) 


(km s" 1 ) 


(deg) 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


11 


NGC 4064 


12 04 11.2 


18 26 36 


SBc(s): 


SB(s)a:pec 


T/C 


12.30 


70 


150 


931 ± 12 


8.8 


NGC 4293 


12 21 12.8 


18 22 57 


Sa pec 


(R)SB(s)0/a 


T/A 


11.20 


67 


66 


930 ± 14 


6.4 


NGC 4351 


12 24 01.6 


12 12 18 


Sc(s) II.3 


SB(rs)ab: pec 


T/N[s] 


13.0 4 


47 


80 


2317 ± 16 


1.7 


NGC 4424 


12 27 11.5 


09 25 15 


Sa pec 


SB(s)a: 


T/C 


12.32 


60 


90 


440 ± 6 


3.1 


NGC 4429 


12 27 26.4 


11 06 29 


S0 3 (6)/Sa pec 


SA(r)0+ 




10.9 


02 


90 


1127 ± 31 


1.5 


NGC 4450 


12 28 29.3 


17 05 07 


Sab pec 


SA(s)ab 


T/A 


10.93 


46 


175 


1958 ± 6 


4.7 


NGC 4569 


12 36 49.8 


13 09 46 


Sab(s) I-II 


SAB(rs)ab 


T/N[s] 


10.2 5 


64 


23 


-232 ± 22 


1.7 


NGC 4580 


12 37 48.6 


05 22 06 


Sc/Sa 


SAB(rs)a pec 


T/N[s] 


12.49 


45 


158 


1036 ± 7 


7.2 


NGC 4606 


12 40 57.6 


11 54 44 


Sa pec 


SB(s)a: 


T/C 


12.69 


67 


40 


1655 ± 16 


2.5 


NGC 4651 


12 43 42.6 


16 23 36 


Sc(r) I-II 


SA(rs)c 


N 


11.36 


51 


71 


804 ± 10 


5.1 


NGC 4694 


12 48 15.1 


10 59 00 


Amorph 


SBO pec 


T/N 


12.19 


42 


140 


1177 ± 11 


4.5 



Note. — (1) Galaxy name; (2) Right Ascension in hours, minutes, and seconds; (3) Declination in degrees, minutes, and seconds; (4) Hubble 
Types from BST, Sandagc & Tammann 1987, or Sandagc Sz Bcdkc 1994; (5) Hubble type from RC3; (6) Star formation class from Koopmann 
& Kcnncy 2004; (7) The total, face-on blue magnitude from RC3; (8) Inclination from Koopmann et al. 2001; (9) Optical P.A from Koopmann 
ct al. 2001; (10) Heliocentric radial velocity from HypcrLEDA; (11) The projected angular distance in degrees of thegalaxy from M87. 
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Table 2 

NIR Properties of Sample Galaxies 







H c 


-R20.5 


Name 


mag 


arcsec 


(1 


) 


(2) 


(3) 


NGC 


4064 


8.78 


64.0 a 


NGC 


4293 


7.35 


156.30 


NGC 


4351 


10.44 


35.50 


NGC 


4424 


9.17 


75.69 a 


NGC 


4429 


6.76 


166.50 


NGC 


4450 


6.90 


147.4 


NGC 


4569 


6.77 


195.00 


NGC 


4580 


8.77 


58.79 


NGC 


4606 


9.30 


60.93 


NGC 


4651 


8.25 


69.06 


NGC 


4694 


9.03 


61.30 



a 2MASS Galaxy Atlas, Jar- 
rctt et al. 2003. 

Note. — (1) Galaxy name; 

(2) Apparent magnitude in ri- 
band (Gavazzi et al. 1999); 

(3) Radius where the surface 
brightness is 20.5 mag arcsec" 2 
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Table 3 

Best Model Parameters of Sample Galaxies 





inc 








aa 


a 








Name 


(°) 


a/ao 




T© Mpc 


T© Mpc 




Comments 




(2) 


(3) 








f5 % l 






(6) 


NGC 4064 


70 


1.82 ± 


0. 


.26 


32 


58 ± 8 


Major 


axis. 


, mean x > 5" 


NGC 4293 


70 


2.68 ± 


0. 


.72 


32 


86 ± 23 


Major 


axis. 


, mean x > 5" 


NGC 4351 


47 


1.12 ± 


0. 


.13 


32 


36 ± 4 


Major 


axis. 


, mean x > 5" 


NGC 4424 


72 


1.17 ± 





33 


32 


37 ± 11 


Major 


axis. 


, mean x > 5" 


NGC 4429 


66 


3.0 ± 





.5 


32 


96 ± 16 


Major 


axis, 


mean x > 12" a 


NGC 4450 


50 


2.02 ± 





.34 


32 


65 ± 1 


Major 


axis. 


, mean x > 5" 


NGC 4569 


66 


1.14 ± 





.09 


40 


46 ± 4 


Major 


axis. 


, mean x > 5" 


NGC 4580 


46 


2.19 ± 





.09 


32 


70 ± 3 


Major 


axis. 


, mean x > 5" 


NGC 4606 


68 


1.35 ± 


0. 


.15 


32 


43 ± 5 


Major 


axis. 


, mean x > 5" 


NGC 4651 


51 


1.53 ± 


0. 


.15 


55 


84 ± 8 


Major 


axis. 


, mean x > 5" 


NGC 4694 


59 


0.43 ± 


0. 


.09 


32 


14 ± 3 


Major 


axis. 


, mean x > 5" 



a Rcgion dominated by circumnuclear disk was avoided 



Note. — (1) Galaxy name; (2) Inclination used for obtaining oblate models; (3) Mean 
a/ao ratio; (4) Model ao] (5) Best a parameter; (6) Comments. 



Table 4 

Magnitudes and Distance Estimates 





M H 


H 


^stellar 


T 


e 




dm 






Name 


mag 


mag 


Mpc 


To 


(°) 


Mpc 


Mpc 


km s 


^stellar/dill 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NGC 4064 


-22.5 ±° l 


3l.3±g;? 


18.0 ±J;| 

14.1 ±;J 


3.2 ±° o[ l 

6-1 ±1:1 


8.8 


2.9 ±°;° 


9.8 ± 0.4 


929 ± 3 


1.84 ±° ™ 


NGC 4293 


-23.4 ±g;l 


30.8 ±°:l 


6.4 


3.2 ±1:1 


14.2 ±J j 


926 ± 4 


0.99 ±8 : S 


NGC 4351 


-21.1 ± 0.1 


31.5 ± 0.1 


20.3 ±°2 


1.8 ± 0.2 


1.7 


o a _1_0.8 

3.6 ±!. 8 


11.2 ±l l 


2310 ± 2 


1-8 ± B oi 


NGC 4424 


-21.7 ± 0.3 


30.9 ± 0.3 


15.2 ± 1.9 


2.4 ± 0.8 


3.1 


1.8 ± 1.6 


4.1 ± 0.2 


442 ± 4 


3.7 ± 0.7 


NGC 4429 


-24.3 ± 0.1 


31.1 ± 0.1 


16.3 ± 0.9 


7.0 ± 1.9 


1.5 


0.7 ± 0.7 




1094 ± 6 




NGC 4450 a 


-23.6 2 


30.5 ±%l 


12.6 ±J;° 


5.9 ± 1.0 


4.7 


4.4 ±0-6 


14.1 ±26 


1953 ± 2 


0.89 ±° .il 


NGC 4569 b 


-23.20 ±°';j£ 


29.97 ±° °l 


9.9 ± 0.2 


4.7 ± 0.4 


1.7 


7.0 ± 0.2 


9.4 ±\i 


-222 ± 6 


1-1 ±o:b 


NGC 4580 


-22.95 ± 0.05 


31.72 ± 0.05 


22.1 ± 0.5 


3.2 ± 0.2 


7.2 


5.8 ± 0.5 


19.3 ± 1.1 


1031 ± 4 


1.1 ± 0.1 


NGC 4606 


-22.20 ±°-°l 


31.50 ±°;Jjj 


19.9 ±i;« 


2.16 ± 0.3 


2.5 


3 25 ± 10 


i o 7 + D - 2 

8-3 


1636 ± 5 


1.6 ± 0.1 


NGC 4651 


-23.90 ±° },1 


32.15 ±g;« 


26.9 ±12 


3.1 ± 0.3 


5.1 


10.3 ±\ l 


22.3 ±2;? 


784 ± 2 


1.2 ± 0.4 


NGC 4694 


-21.6 ± 0.2 


30.6 ± 0.2 


13.4 ±i 3 


1.1 ± 0.3 


4.5 


3.6 ±° | 




1174 ± 22 





a If wc add a dark matter halo in order to match the Ha rotation curve we have that Mh — -23.8 ± 0.1 mag, fi — 30.7 ± O.lmag , and 
dsteiiar - 13.8 ±g;| Mpc 

b If we add a dark matter halo in order to match the Ha rotation curve we have that Mh — -23.8 ± 0.01 mag, fi — 30.57 ± 0.01 mag , 
and ^stellar = 13.0 ± 0.1 Mpc 

Note. — (1) Galaxy name; (2) Absolute magnitude in H-band; (3) Distance modulus; (4) SKB distance; (5) Mass-luminosity ratio in the 
R-band; (6) Angular distance to M 87; (7) Barycentric SKB distance to M 87; (8) HI based distance; (9) Stellar systemic velocity; (10) 
Ratio between dc SKB distance and HI based distance. 
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